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Abstract 

We will discuss how the decays of charm mesons can be used to study 
light mesons spectroscopy, by presenting recent results of Dalitz plot anal- 
yses using data from Fermilab experiment E791. Emphasis will be on 
scalar mesons, which are found to have large contribution to the D decays 
studied. In addition to the usual extraction of decay fractions and relative 
phases of the intermediate amplitudes, the Dalitz plot technique is used 
to measure masses and widths of scalar resonances. From the D s decay, 
we obtain masses and widths of /o(980) and /o(1370). We find evidence 
for a light and broad scalar resonance, the a meson, in D + — > tv~tt + tv + 
decay. Preliminary studies also show evidence for a light and broad reso- 
nance, the k meson, in D + — > K~tx + -k + decay. These results illustrate the 
potential of charm decays as a laboratory for the study of light mesons. 

1 Introduction 

The decays of D mesons can be seen as a new environment for the study of light 
meson physics, due to the high quality of the current samples, the fact that the 
initial state is always well defined (defined mass, CP state), and usually with 
a small non-resonant contribution. Particularly, these decays can shed light on 
the scalar mesons, some of them long-standing sources of uncertainty. 

Here we present results for the Dalitz plot analyses of the decays Df — ► 
7r~7r + 7r + HJ and D + — > ir~ir + ir + and also preliminary results for the decay 
D + — ► K~ir + ir + , using data from the Fermilab E791 experiment]^]. 

In the Dj — > tt~t: + it + decay, previous measurements showed the dominance 
of scalar resonant states^, |). F° r instance, the decay /o(980)7r + seems to 
contribute most. Nevertheless, the nature of the /o(980) is still a puzzle, and in 
particular its width is poorly measured|(|. From the E791 Dalitz plot analysis 

1 E791 ran in 91-92 with a 500 GeV/c tt~ beam. See 



of the D+ — + 7r _ 7r + 7T + decays, we are able to obtain measurements for the mass 
and width of this state, as well as for the /o(1370). For the D + — * ir~ir + ir + 
decay, we are not able to get a good description of the data using the same 
intermediate states as for — > ir~ir + ir + . We find strong evidence for a new 
light and broad scalar resonance, the <r(500), and measure its mass and width. 
The D + —* o"7r + channel accounts for about half of the total D + — > h~ii + ti + 
decay rate. Other experiments have presented inconsistent evidence for low- 
mass 7T7T resonances in partial wave analyses @ , resulting in ambiguous results 
for the characteristics of such particles || . 

For the D + — > K~ir + ir + decay, we present preliminary results from a high 
statistics sample (almost 23,000 decays). Like previous results for this channel || 
[Tof we find that the non-resonant decay is dominant, which is unusual in D 
decays, but also find important discrepancies between fit model and data. By 
including a new scalar state, with unconstrained mass and width, we get a 
much better description of the data. This resonance appears as light and broad 
according to the fit. Present discussions about the existence of such a state, 
referred to as the n in the literature, are very controversial Q. 

2 D + ,D+ —> 7T 7r + 7r + Dalitz plot analyses 

In Fig. |l| we show the tt^tt + tt + invariant mass distribution for the sample col- 
lected by E791 after reconstruction and selection criteria || . Besides combina- 
torial background, reflections from the decays D + — > K ~tt + it + , D + — > K~tt + 
(plus one extra track) and Df — » rj'ir + , r} 1 — * p° (770)7 are taken into ac- 
count. The hatched regions in Fig. |l| show the samples used for the Dalitz 
analyses. There are 937 and 1686 candidate events for Df and D + respectively, 
with a signal to background ratio of about 2:1. 

The Dalitz plot corresponding to Df — > ir~Tr + ir + events is shown in 
Fig. |](a). The narrow horizontal and vertical bands at S12 = m 2 (7rj~7r^) and 
Si3 = m 2 (7rf 7r^") just below 1 GeV 2 /c 4 correspond to the /o(980)7r + state. At 
the upper edge of the diagonal, there is another concentration of events centered 
at S12 — Si3 ~1.8 GeV 2 /c 4 , corresponding to the /2(1270)7r + , / (1370)7r + , and 
/3°(1450)7r + contributions. From the D + — > 7r _ 7r + 7r + Dalitz plot (Fig. §(b)) we 
see clearly the bands corresponding to the p°(770)7r + and /o(980)7r + channels, 
and an excess of events at low tt~it + invariant-mass squared. 

To study the resonant structure of these decays, unbinncd maximum like- 
lihood fits are applied to the samples. The Dalitz plot distributions are fitted 
to a sum of signal and background probability distribution functions (PDFs). 
For each candidate event, the signal PDF is written as the square of the to- 
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Figure 1: The 7r~7r + 7r + invariant mass spectrum. The dotted line represents 
the D° — > A" _ 7r + plus Df — > ?/7r + reflections and the dashed line is the total 
background. Events used for the Dalitz analyses are in the hatched areas. 
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Figure 2: (a) The Df — > 7r 7T+7T+ Dalitz plot and (b) </ie Z) + 7r 7r+7r^ 
Dalitz plot. Since there are two identical particles, the plots are symmetrized. 



tal physical amplitude, which is constructed from a coherent sum of Lorentz 
invariant amplitudes corresponding to the various resonant channels plus the 
non-resonant decay. The signal PDF is weighted by the acceptance across the 
Dalitz plot (obtained by Monte Carlo (MC)) and by the level of signal to back- 
ground for each event, as given by the line shape of Fig. [|. 

The resonant channels we include in the fits are p°(770)7r + , /o(980)7r + , 
/ 2 (1270)7r+, /o(1370)7r+, and p (1450)tt+. We assume the non-resonant am- 
plitude to be uniform across the Dalitz plot. Each resonant amplitude, except 
that for the /o(980), is parameterized as a product of form factors, a relativis- 
tic Breit-Wigner (B-W) function, and an angular momentum amplitude which 
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depends on the spin of the resonance, 

, _ Fp Fr M ( n ,] r(m lz )-r u m ° (p*\ 2J+1 f Up*) (1) 

mf 2 - ml + im r(TOi 2 ) ' 12 ° m 12 \pqJ Fg(Po) 

Above, 77Ji2 is the 7r + 7r~ invariant mass for the candidate resonance. Since 
there are two like-charge pions, each signal amplitude is Bose-symmetrized, 
A n = Ai [(12)3] +Ai [(13)2]. The quantities F D and F R are the Blatt-Weisskopf 
damping factors [O respectively for the D and the resonancen, p* is the pion 
momentum in the resonance rest frame at mass 777,12 (Po = P*("7o))- -A4n 
describes the angular distribution due to the spin J of the resonance. 

For the /o(980)7r + we use a coupled-channel B-W function, following the 
parameterization of the WA76 Collaboration|13| , 



1 

-ml + imo(T n + T K ) 



BW M 980) = —2 _ 2 I ZZ rF , t^__n ' ( 2 ) 



r n = g^y/m^/A-ml, T K = ^- (^J m 2 ^ / A ~ m 2 K+ + yj / '4 - m 2 K ^j . 

(3) 

We multiply each amplitude by a complex coefficient, c n = ane 5 ". The fit 
parameters are the magnitudes, a n , and the phases, S n , which accommodate the 
final state interactions. 



2.1 Results for Df 
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The Df — > 7r~7r + 7r + Dalitz plot (Fig. 0(a)) is fitted to obtain not only the 
relative contributions (and phases) of the possible sub-channels, but also the 
parameters of the /o(980) state, g^, gx, and too, as well as the mass and width 
of the /o(1370). That is, they are determined directly from the data, as free 
parameters in the fit. The other resonance masses and widths are taken from 
the PDG[^]. The resulting fractions and phases are shown in Table [l]. 

The measured /o(980) parameters are m = 977 ± 3 ± 2 MeV/c 2 , g* = 0.09 
± 0.01 ± 0.01 and g K = 0.02 ± 0.04 ± 0.03. Our value for g^ is in very good 
agreement with OPAL and MARKII results Q, but WA76 Q found a much 
larger value, g n — 0.28 ± 0.04. Our value of gx indicates a small coupling 
of /o(980) to KK. The values of the Jo (980) mass and of g n , as well as the 
magnitudes and phases of the resonant amplitudes, are relatively insensitive to 
the value of gK- Both OPAL and MARKII results are also insensitive to the 
value of gK- WA76, on the contrary, measured gK = 0.56 ± 0.18. 

2 For three pions analyses, the effective radii for and Fr are set to 3.0GeV — 1 ; for 
the D+ — » i<r _ 7r+7r+ analysis of section H, the radii taken for Fn and Fr are respectively 
S.OGeV" 1 and l.SGeV" 1 . 
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Table 1: Dalitz fit results for — > n n + 7r' 



Decay Mode 


Phasc(°) 


Fraction(%) 


/o(980)tt+ 
non-rcson. 
p°(770)7r+ 
/ 2 (1270)tt+ 
/o(1370)7r+ 
P (1450)tt + 


(fixed) 
181 ±94 ±51 
109 ± 24 ± 5 
133 ± 13 ±28 
198 ± 19 ±27 
162 ±26 ±17 


56.5 ±4.3 ±4.7 
0.5 ± 1.4 ± 1.7 
5.8 ±2.3 ±3.7 
19.7 ±3.3 ±0.6 
32.4 ± 7.7 ± 1.9 
4.4 ±2.1 ±0.2 
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Figure 3: S12 and S13 (m^) projections for — > ir ir + ir + data (dots) and 
our best fit (solid). The hashed area corresponds to background. 

We have also fit the Dalitz plot using for the /q(980) the same B-W function 
as for the other resonances. We find m = 975 ± 3 MeV/c 2 and T = 44 ± 2 ± 2 
MeV/c 2 , and the results for fractions and phases are indistinguishable. 

A x 2 distribution is produced from the difference in densities for model 
(from a fast-MC algorithm) and data. From the x 2 an d the number of degrees 
of freedom (z/), the confidence level of the model is 35% |IJ. A visual comparison 
between the fit model and the data can be seen in Fig. [3] where we show the 
sum of the S12 and S13 projections for data (points) and model (solid lines, from 
fast-MC). 

If the Df — ► 7r~7r + 7r + decay is dominated by the Cabibbo-favored spectator 
mechanism, we would expect final states with a large ss content. As we can 
see by the results of Table [j], approximately half of the Df — > tt^tt + tt + rate is 
produced via /o(980)7r + . If the spectator amplitude is dominant in this decay, 
this would support the interpretation of the /o(980) as an ss state. On the other 
hand, the large contribution from the intermediate state /o(1370)7r + indicates 
the presence of either VF-annihilation amplitudes or strong rescattering in the 
final state. In fact, the /o(1370)7r + is not observed in the — > K + K~ir + 
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final state [fl5|, pointing to the /o(1370) being a non-ss particle, as suggested by 
the naive quark modelQ. There is no evidence in the decay for a low-mass 
broad scalar particle as seen in the D + decay, discussed below. 

2.2 Results for D + -k--k + -k + 

In a first approach, we try to fit the D + — > it~-K + it + Dalitz plot (Fig. |^(b)) using 
the same amplitudes used for the Df — > n~Tr + ir + analysis. We find that, with 
this model, the non-resonant, the p°(1450)7r + , and the p°(770)ir + amplitudes 
dominate (results shown on the first column of Table ^|), reproducing the same 
qualitative features from those reported previously j|, However, this model 
does not describe the data satisfactorily, especially at low tt~~it + mass squared, 
as can be seen from Fig. ||(a). The x 2 / v obtained from the binned Dalitz plot 
for this model is 1.6, with a confidence level less than 1CP 5 . 

To investigate the possibility that another it~it + resonance contributes to 
the D + — > 7r~7r + 7r + decay, we add an extra scalar resonance amplitude to the 
signal PDF. We allow its mass and width to float as free fit parameters. 

We find that this model improves our fit substantially, converging to values 
of mass and width of this scalar resonance, called ct(500), of 478±23±17 MeV/c 2 
and 324+ 4 q ±21 MeV/c 2 , respectively. According to this model, this amplitude 
produces the largest decay fraction, as shown in the second column of Table 
||; the non- resonant amplitude, which is dominant in the model without air + , 
drops substantially. This model describes the data much better, as can be seen 
by the irir mass squared projection in Fig. ||(b). The x 2 / v is now 0.9, with a 
corresponding confidence level of 91%. 

To better understand our data, we also fit it with vector, tensor, and toy 
models for this extra amplitude, allowing the masses, widths, and relative ampli- 
tudes to float freely. The vector and tensor models test the angular distribution 
of the signal. The toy model tests the phase variation expected from a B-W 
amplitude by forcing this amplitude to have a constant relative phase but still 
a B-W shape vs. mass. All these alternative models fail to describe the data as 
well as the scalar (regular) B-W amplitude (see Q). 

3 Preliminary results for D + — > K ~7r + 7r + Dalitz 
analysis 

Previous analyses of the Dalitz plot of the decay D + — » i^ _ 7r + 7r + |j| [l(| showed 
the unusual dominance of the non-resonant decay. The most recent result comes 
from E687 [^0) with a sample of about 8000 events, where besides the non- 



G 



Table 2: Dalitz fit results for D + — > tt tt + tt + . First errors are statistical, 
second systematics (only for fit with er7r + mode). 



Decay 


Fit without o"7T + 


Fit with f77r + 


Mode 


Phase (°) 


Fraction(%) 


Phase(°) 


Fraction(%) 


(77T + 






206 ± 8 ± 5 


46.3 ±9.0 ±2.1 


P°(770)7T+ 


(fixed) 


20.8 ±2.4 


(fixed) 


33.6 ±3.2 ±2.2 


non-reson. 


150 ± 12 


38.6 ±9.7 


57 ± 20 ± 6 


7.8 ±6.0 ±2.7 


fo (980)^+ 


152 ± 16 


7.4 ± 1.4 


165 ±11 ±3 


6.2 ± 1.3 ±0.4 


/ 2 (1270)7r+ 


103 ± 16 


6.3 ± 1.9 


57±8±3 


19.4 ±2.5 ±0.4 


/o(1370)7r+ 


143 ± 10 


10.7 ± 3.1 


105±18±1 


2.3 ± 1.5 ±0.8 


p°(1450)7r+ 


46 ± 15 


22.6 ±3.7 


319 ±39 ±11 


0.7 ±0.7 ±0.3 




(a) (b) 

Figure 4: s\2 and S13 (m~ w ) projections for D + — > data (dots) and 

our best fit (solid) for models (a) without and (b) with o~n + amplitude. The 
dashed distribution corresponds to the expected background level. 

resonant (about 95% of the decay rate) they find contributions from the channels 
i^(1430)7r+, jf*(890)7r+, and K*(1680)tt+. However, they don't obtain a good 
description of the data, arriving to a x 2 jv of 3. 

In Fig. |E|(a) we show the K~ir + n + invariant mass distribution for the sample 
collected by E791 after reconstruction and selection criteria. Besides combinato- 
rial background, the other main source of background comes from the reflection 
of the decay D+ -> K-R+-K+ (through K*K+ and (jm + ). The total level of 
background is shown by the filled area in Fig. |(a). The hashed area corre- 
sponds to the sample used for the Dalitz analysis. There are 22890 events in 
this sample, where about 6% correspond to background. 

The Dalitz plot of the signal-region events is shown in Fig. ||(b). The plot 
presents a rich structure, where we can observe the clear bands from K* (890)ir + , 
and an accumulation of events at the upper edge of the diagonal, due to the 
heavier resonances. 
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Figure 5: (a) The Kmr invariant mass spectrum. The shaded area is back- 
ground; (b) Dalitz plot corresponding to the events in the dashed area of (a). 



The formalism used for the D + — * K~ir + ir + Dalitz plot fit is essentially the 
same as for the D + , Df three pions decays^. The results presented here are all 
preliminary, in particular no studies of systematic effects are addressed. 

All possible known Kn resonances [[| are included in the model in a first step. 
We find significant contributions from the same channels observed previously || 
|Io| , but with our higher statistics sample we also measure a small but significant 
contribution from ^|(1430)7r + . The decay fractions and phases for this model 
are shown in Table [|, first column. We confirm the high non-resonant rate 
(over 100%); there is an impressive interference pattern according to this model, 
reflected by the sum of the fractions being around 150%. 

However, this model does not provide a good description of the data. The 
calculated x 2 / l/ horn the binned Dalitz plot distribution is 2.7. By observing 
the m 2 (Kiti ow ) and m 2 (KTthigh) projections for data (points) and model (solid) 
in Fig. |6|(a), we see discrepancies at low Kir mass squared (below 0.6 GeV 2 /c 4 , 
where data has very small error bars) and also near 2.5 GeV 2 /c 4 . 

The strong evidence for the cr(500) from the D + — > 7r~7r + 7r + decay leads 
us to question whether a new resonance state could be responsible for the dif- 
ficulty in modeling the D + — > K~n + ir + data. The light and broad scalar Ktt 
resonance, the k (a c-meson nonet member), is presently the subject of many 
discussions, with controversial results from Kn scattering analyses [pi[ [PI) . 

3 Besides the differences in the factor Fr addressed before, here we define the relativistic 
Breit-Wigner with a factor of (-1) with respect to equation ij for easier comparison to previous 



analysis ]1C 
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Table 3: Dalitz fit results for D+ -> R--k + t: + . PRELIMINARY; statistical 
errors only. 



Decay 


Fit without kit + 


Fit with ktt + 


Mode 


Phase(°) 


Fraction(%) 


Phase(°) 


Fraction(%) 


non-reson. 


(fixed) 


103.9 ±2.3 


(fixed) 


52.4 ±8.5 








187 ± 11 


20.7 ±5.4 


^*(890)tt+ 


51 ±1 


12.7 ±0.4 


15 ±6 


11.3 ±0.4 


^(1430)tt+ 


63 ±1 


33.6 ± 1.4 


63 ±5 


17.6 ± 1.7 


if 2 *(1430)7r+ 


51 ±6 


0.5 ±0.1 


348 ± 10 


0.3 ±0.1 


A*(1680)tt+ 


70 ±3 


3.8 ±0.3 


53 ±6 


3.3 ±0.7 



We introduce in our model an extra spin-0 resonant amplitude. Both mass 
and width of this hypothetical k state are allowed to float in the fitting pro- 
cedure. Moreover, since the mass and width of the other scalar resonance, 
Aq(1430), were obtained from the LASS experiment in the absence of this ex- 
tra amplitude fl6|, we also allow these parameters to float. The fit converges 
to values of the k mass and width of 815 ± 30 MeV/c 2 and 560 ± 116 MeV/c 2 
respectively, and to values for the Aq(1430) mass and width of 1465 ±6 MeV/c 2 
and 182 ± 10 MeV/c 2 respectively. The results for k are consistent to this state 
being a light and broad resonance [jll). The results for the Aq(1430) show this 
state being narrower and heavier than obtained by LASS jfj, [l(J. The decay 
fractions and phases from this model are shown in Table |^, second column. The 
non-resonant fraction drops to about 50% and the ktt + accounts for about 20% 
of the decay, representing the second main contribution, followed by the other 
scalar state, Kq (1430)7r + . The fit quality of this model is considerably superior 
to the model without the ktt + state: the x 2 l v is now 0-9, the confidence level 
of the fit is 79%. As can be seen by the projections in Fig. |6|(b), there is very 
good agreement between model and data. 

As was done for testing the a state, we test a variety of other models to 
check their ability to explain the data. We represent the extra amplitude by 
vector, tensor and toy (B-W with no phase variation) models, allowing mass 
and width to float. None of these models reproduced the data as well as the 
scalar resonance hypothesis. 

4 Conclusions 

Dalitz plot analyses presented here from E791 data point to the importance of 
the scalar resonant states in the 37T and Ktttt final states of D mesons. From 
the Df — > 7r~7r + 7r + analysis, we find that about 90% of the decay rate is 
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due to /o7r + intermediate states. Almost half of the D + — > h~it + ti + decay 
rate proceeds through a new scalar state tr(500)7r + . From D + — > K~n + it + 
decay, we present evidence for a new scalar state ktt + and, although about 50% 
of the decay rate corresponds to the non-resonant channel, 40% comes from 
the scalar states. We are able to measure the masses and widths of these scalar 
resonances. These results illustrate the potential of many-body D meson decays 
for the study of light-meson spectroscopy. 



References 

[1] E.M. Aitala etal. (E791 Coll.), |hep-ex/0007027| , accepted by PRL. 
[2] E.M. Aitala etal. (E791 Coll.), |hep-ex/0007028| , accepted by PRL. 
[3] J.Appel, Ann. Rev. Nucl. Part. Sci. 42, 367 (1992); D.Summers etal., 



hep-ex/0009015; S.Amato etal., Nucl. Instr. Meth. A 324, 535 (1993); 



E.M. Aitala etal. (E791 Coll.), Eur. Phys. J. direct C 4, 1 (1999). 
[4] J.C.Anjos etal. (E691 Coll.), Phys. Rev. Lett. 62, 125 (1989). 
[5] P.L.Frabetti etal. (E687 Coll.), Phys. Lett. B 407, 79 (1997). 



10 



[6] C.Caso etal. (Particle Data Group), Eur.Phys. J.C 15, 1 (2000). 

[7] D. Barberis etal. (WA102 Coll.), Phys. Lett. B 453, 316 (1999); D.M. Asner 
etal. (CLEO Coll.), Phys. Rev. D 61, 012002 (2000); D. Aide, etal. (GAMS 
Coll.), Phys. Lett. B 397 , 350 (1997). 



N. Tornqvist, |hcp-ph/9904346j . 



[9] J.C.Anjos etal. (E691 Coll.), Phys. Rev. D 48, 56 (1993). 
[10] P.L.Frabetti etal. (E687 Coll.), Phys. Lett. B 331, 217 (1994). 
[11] D. Black etal, Phys. Rev. D 58, 054012 (1998); S. Ishida etal, Prog. 



Theor. Phys. 98, 621 (1997); S.N. Cherry and M.R. Pennington, frep- 
ph/0005208 ; and references therein. 



[12] J.Blatt and V.Weisskopf, Theoretical Nuclear Physics (Wiley, 1962). 

[13] T.A.Armstrong etal. (WA76 Coll.), Z. Phys. C 51, 351 (1991). 

[14] K.Ackerstaff etal. (OPAL Coll.), Eur. Phys. J.C 4, 19 (1998); G. Gidal 
etal (MARKII Coll.), Phys. Lett. B 107, 153 (1981). 

[15] P.L.Frabetti etal (E687 Coll.), Phys. Lett. B 351, 591 (1995). 

[16] D. Aston etal. (LASS Coll.), Nucl.Phys.B 296, 493 (1988). 



11 



